Nuclear matrix elements of 0v/3f3 decay with improved short-range correlations 
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Nuclear matrix elements of the neutrinoless double beta (Ov/3/3) decays of 96 Zr, 100 Mo, 116 Cd, 
128 Te, 130 Te and 136 Xe are calculated for the light-neutrino exchange mechanism by using the proton- 
neutron quasiparticle random-phase approximation (pnQRPA) with a realistic nucleon-nucleon 
force. The g pp parameter of the pnQRPA is fixed by the data on the two-neutrino double beta 
decays and single beta decays. The finite size of a nucleon, the higher-order terms of nucleonic weak 
currents, and the nucleon-nucleon short-range correlations (s.r.c) are taken into account. The s.r.c. 
are computed by the traditional Jastrow method and by the more advanced unitary correlation op- 
erator method (UCOM). Comparison of the results obtained by the two methods is carried out. The 
UCOM computed matrix elements turn out to be considerably larger than the Jastrow computed 
ones. This result is important for the assessment of the neutrino-mass sensitivity of the present and 
future double beta experiments. 
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PACS numbers: 21.60. Cs, 23.40. He, 27.60. +j 
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I. INTRODUCTION 

The neutrinoless double beta (Ov/3/3) decay plays a key role in the search for massive Majorana neutrinos and 
their mass scale. The experimental search for 0vf3(3 decay has become front-line physics due to the verification 
of the existence of neutrino mass by the oscillation experiments [H and the claimed discovery of the 0^/5/3 deca- 
ff, la]. At present there are two important experiments taking data, namely the NEMO 3 [4J and CUORICINO 
experiments. A host of important future experiments are under R&D planning and construction @|. For all these 
expensive experiments the computed values of the involved nuclear matrix elements have become an important issue 
0,1, ©, 03 ■ They are essential when one starts to extract quantitative neutrino properties from the measured data. 

Many nuclear models of different types have been deviced to compute the nuclear matrix elements of the Ov/3/3 
dec ay II ill. There are essentially two complementary families of nuclear models on the market, the nuclear shell model 
[T^.ll3l.ll4| and the proton-neutron quasiparticle random-phase approximation (pnQRPA) [TJ, EH, EH] ■ The pnQRPA 
is constructed to describe the energy levels of odd-odd nuclei and their beta decays to the neighboring even-even 
nuclei pj}- Also its derivative, renormalized pnQRPA [l8j], has been used [Isl l20| to compute double-beta matrix 
elements. The pnQRPA (and the renormalized pnQRPA) calculations can be fine tuned by the so-called particle- 
particle strength parameter, g pp , that controls the magnitude of the proton-neutron two-body interaction for the 1 + 
intermediate states in double beta decay [l3, lla]- There are two ways to fix the value of this parameter, either by 
using the data on two-neutrino double beta (2v(3/3) decays [20] or the data on single beta decays [2~ll. 1221] . In this work 
we use the 2v(3(3 data to fix the possible values of g pp and cross check it by single beta decays wherever possible. 
In the mass mode of the Ov/3/3 decay a light virtual Majorana neutrino is exchanged by the two decaying neutrons 
t^. \ of the initial nucleus. The average exchanged momentum is large so that the two neutrons tend to overlap. To 
prevent this a Jastrow type of correlation function was introduced in [l2l . [23| following the parametrization of Miller 
and Spencer [Hj]. This method, although microscopically inspired, is just a phenomenological way to introduce short- 
range correlations into the two-nucleon relative wave function. The Jastrow function simply cuts off the wave function 
of the two nucleons at short relative distances r leading to a violation of the norm of the wave function. 

In the present calculations we improve on the Jastrow method by engaging the more sophisticated microscopic 
approach of unitary correlation operator method (UCOM) [25| . In the UCOM one obtains the correlated many- 
particle state from the uncorrelated one by a unitary transformation and thus the norm of the correlated state is 
conserved and no amplitude is lost in the relative wave function. In the 0v(3f3 calculations this leads to a more 
complete description of the relative wave function for small r, as was demonstrated in [2(| for the decays of 48 Ca and 
76 Ge. In this work and in [13], where the decays of 76 Ge and 82 Se were analyzed, it is demonstrated that the Jastrow 
procedure leads to the excessive reduction of 25% - 40% in the magnitudes of the Ov/3/3 nuclear matrix elements. At 
the same time the UCOM reduces the magnitudes of the matrix elements only by 4% - 16%. The magnitude of the 
short-range corrections affects the magnitudes of the nuclear matrix elements which, in turn, dictate the neutrino-mass 
sensitivity of the potentially succesful future double beta experiments. The notable differences between the Jastrow 
and UCOM corrections influence the cost estimates of large-scale experiments if a given neutrino-mass sensitivity is 
wanted. 

In the article we continue the work of where the Ov/3/3 nuclear matrix elements of 76 Ge and 82 Se were derived. 



We apply the UCOM and Jastrow short-range correlations on matrix elements derived by the pnQRPA method and 
corrected for the higher-order terms of nucleonic weak current and the nucleoli's finite size using the recipes of [2oLl2gj|. 
We analyze the 0i//3/?-decay matrix elements of 96 Zr, 100 Mo, 116 Cd, 128 Te, 130 Te and 136 Xe for all the mentioned 
corrections. The necessary theoretical background is briefly described in Sec. II and the numerical application is 
reviewed in Sec. III. The results are discussed in Sec. IV, and the summary and conclusions are presented in Sec. V. 



II. THEORETICAL BACKGROUND 



We start this short review of the theory by presenting the expression for the half-life of the 2v(3(3 decay: 
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4/ 2 W-0+)j =GM M™, . (1) 

The transition proceeds from the initial ground state, 0+, to the final ground state, 0+. Here G^ 2 ^ is an integral over 

the phase t 
written as 



the phase space of the leptonic variables [Hi, |29(. The involved double Gamow-Teller matrix element, M^q T , can be 



M (2„) = v (QtwE^tmw^) 

x(l+||5>(7)*7ll°i + )' (2) 

3 

where the transition operators are the usual Gamow-Teller operators for 0~ transitions, Qpp is the 2vf}(3 Q value, E n 
is the energy of the nth intermediate state, M\ is the mass energy of the initial nucleus, and m e is the rest-mass of the 
electron. It has to be noted here that the expression ((2J is scaled by the electron rest mass to yield a dimensionless 
matrix element. This definition deviates from that of some other authors, like [23], where the scaling is not done. 

The 0v(3f3 decay can proceed via the exchange of a light virtual Majorana neutrino. Assuming this neutrino-mass 
mechanism to be the dominant one, the inverse of the 0v/3f3 half-life can be written as 

^(O^O+f^G^f^j^M^) 2 , M(°">=Jl4¥- M^+M^. (3) 

In the above expression M^ ^ is the total nuclear matrix element consisting of the Fermi, Gamow-Teller and tensor 
contributions. The effective mass of the neutrino is given by 



(m„) = X>f mj\U ej \ 2 , (4) 

3 



where A^ p is the CP phase and U e j is a component of the neutrino mixing matrix. The definition of the leptonic 

phase-space factor can be found in fill ]. 

The nuclear matrix elements involved in the mass mode of the 0^/3/3 decay are defined as 

= £(0+||Mr WWSJ E a )||0+) ) (5) 

a 

m gt = J2(0f\\h G T(r mn ,E a )a m -(T n \\0f), (6) 

a 

where the summation runs over all the intermediate states and the integration is taken over the relative coordinate 
Tmn = l r m — r„| between the nucleons m and n. The neutrino potential h.K(r mn ,E a ), K = F, GT, is defined as 

h K (r m n,E a ) = -R A [ dq | F ; /0 jo(qr mn ) , (7) 

7T J q + E a - {Ei+ Ef)/2 



where Ra = 1.2^4 1//3 fm is the nuclear radius and jo is the spherical Bessel function. The term h,K{q 2 ) in (7|) includes 
the contributions arising from the induced currents and the finite nucleon size [13, HE|- 



Next we write the nuclear matrix elements explicitly in the pnQRPA framework. They are given by 



J 17 ,k2,J' pp'nn' 



• P j,}(pp':J'\\0 K \\nn':J') (8) 



Jn' Jp 
.t 7, 



x(o f + || [£^J J^XJ&KK^II L4 5 «JJI°i ). 

where fci and fc 2 label the different pnQRPA solutions for a given multipole J w , The operators Ok inside the 
two-particle matrix element derive from J5]) and and they can be written as 

F = h F (r,E k ), O gt = h GT (r,E k )eT 1 -(T 2 , r = \r x - r 2 | . (9) 

The expression for the pnQRPA transition densities (0^||[Cp,c„']j|| J k± ) and (Jfc 2 ||[c£c„]j||0+), and for the overlap 
factor (JfcJJ^,) can be found e.g. in (ill. \vfy. 

The traditional way to include the short-range correlations in the Ois/3/3 decay calculations is by introducing the 
Jastrow correlation function /j (r) , which depends on the relative distance r = |ri — r 2 | of two nucleons. In the Jastrow 
scheme the uncorrelated operator O is replaced by the correlated operator Cj by a procedure 

(0+||O||0+) - (O+H0JHO+) = (0;||/jO/j||0+) . (10) 

A typical choice for the function /j in 0vf3f3 calculations is (3(J 

/j(r) = l-e- ar (l-6r 2 ) , (11) 

with a = 1.1 fm and b = 0.68 fm 2 . Such application of the Jastrow correlation function is very rudimentary and, 
consequently, the Jastrow correlation cuts out the r < 1 fm part from the relative two-particle wave function. In Ref. 
[2^ | it was demonstrated that this leads to overestimation of the effects of short-range correlations on the many-body 
wave function. 

To circumvent the difficulties associated to the Jastrow correlations the more refined unitary correlation operator 
method (UCOM) [25[ was used in the Qvj3(3 decay calculations of Ref. In UCOM one obtains the correlated 

many-body state from the uncorrelated one as 

l*>=C|*>, (12) 

where C is the unitary correlation operator. The operator C is a product of two unitary operators: C — CqC t , where 
Cfi describes short-range tensor correlations and C r central correlations. Due to the unitary of the operator C the 
norm of the correlated state is conserved. Moreover, since the correlated matrix element of the operator O can be 
written as 

<*|0|*') = (*|C t e>C|5'') = (tf (13) 

it is therefore equivalent to use either correlated states or correlated operators. For the Fermi and Gamow- Teller 
0v(3fi nuclear matrix elements the effect of the tensor correlation operator Cq vanishes and one is thus left only with 
the central correlation operator. The UCOM parameters used in our 0^/3/3 calculations are the Bonn-A parameters 
taken from |3J]|- 



III. NUMERICAL APPLICATION 

In [13] we applied the pnQRPA to compute the 0v(3f3 nuclear matrix elements of 76 Ge and 82 Se in the model space 
containing the lp-0f-2s-ld-0g-0hn/2 single-particle orbitals, both for protons and neutrons. Here we add to this model 
space the spin-orbit partner 0hg/ 2 of the 0h n / 2 orbital, both for protons and neutrons, to describe the decays of 96 Zr 
and 100 Mo. For the rest of the decays, namely for the 116 Cd, 128 Te, 130 Te and 136 Xe decays we have extended the 
proton and neutron model spaces to include the lp-0f-2s-ld-0g-2p-lf-0h single-particle orbitals. The single-particle 
energies were obtained from a spherical Coulomb-corrected Woods-Saxon potential with a standard parametrization, 
optimized for nuclei near the line of beta stability. Slight adjustments were done for some of the energies at the vicinity 



TABLE I: Values of the g pp parameter extracted from the data. First column shows the decay and the second column the 
matrix element values extracted from the 2i//3/3-decay data by [20]. Third column gives the range of g pp corresponding to the 



matrix elements of the second column. 


The g pp ranges 


of the last column were extracted from the available (5 


decay data. 


Decay 


m.e. (2^/3/3) 






ao Zr -» yo Mo 


0.026 - 0.112 


1.06 - 1.11 


- 


100 Mo -» 100 Ru 


0.107-0.181 


1.07 - 1.09 


1.07-1.08 


116 Cd^ 116 Sn 


0.058 - 0.102 


0.97-1.01 


0.82 - 0.84 


128 Te — ► 128 Xe 


0.011-0.037 


0.89 - 0.92 


0.86 - 0.88 


130 Te 130 Xe 


0.014 - 0.054 


0.84 - 0.90 




136 Xe ^ 136 Ba 


< 0.023 


> 0.74 





of the proton and neutron Fermi surfaces to reproduce better the low-energy spectra of the neighboring odd- A nuclei 
and those of the intermediate nuclei. 

We have used the Bonn-A G-matrix as the two-body interaction and we have renormalized it in the standard 
way [H, by fitting the pairing parameters of the BCS by comparing with the phenomenological pairing gaps, 
extracted from the atomic mass tables. The particle-hole parameter g p ^ of the pnQRPA affects the position of 
the giant Gamow- Teller resonance and its value was fixed by the available data on the location of the giant state. 
Due to this phenomenological renormalization of the two-body interaction we did not perform an additional UCOM 
renormalization (34j . 

After fixing all the Hamiltonian parameters the only free parameter left is the proton-neutron particle-particle pa- 
rameter g pp of the pnQRPA. We obtained the physical values of <7 PP by using the method of [2fjLl27j|. Consequently, we 
used the extracted experimental matrix elements of [2(J that include the experimental error limits and the uncertainty 
in the value of the axial-vector coupling constant 1.0 < #a < 1.25. The resulting intervals of the experimental 2v/3{3 
matrix elements are shown in the second column of Table HI and they are scaled by the electron rest mass according to 
([2]). By performing the pnQRPA calculations of the 2v(3(3 matrix elements the ranges of experimental matrix elements 
were subsequently converted to the intervals of g pp values shown in column three of Table [B For some cases there 
exists log ft data on (3~ decay from the first 1 + state of the intermediate nucleus to the ground state of the double 
beta daughter nucleus. Applying the above procedure to this data leads to the experimental j3~ matrix elements and 
the corresponding ranges of <? pp , listed in the last column of Table |H 

As can be seen from Table U the values of g pp extracted from the (5~ and 2v(i(i data are (roughly) compatible for 
the decays of 100 Mo and 128 Te, whereas for the 116 Cd decay this is not the case. This discrepancy was already pointed 
out in [2l|. The reason for this discrepancy is not clear but an interesting observation is that the 116 Cd decay obeys 
very closely the single-state dominance hypothesis whereas the 128 Te decay obeys it slightly less and the 100 Mo decay 
the least, as clearly shown in Table 1 of [351 ] . 



IV. DISCUSSION OF THE RESULTS 

In Table HI we show the evolution of the values of the 0v(if3 nuclear matrix elements as we add more corrections to 
the bare matrix element. In the first column we list the used g pp value which was taken to be in the middle of the 
g pp interval of Table U The following columns list the bare matrix element (b.m.e.), the matrix element including 
the higher-order terms of the nucleonic weak current (b.m.e.+A), and the matrix element with finite nucleon size 
effects added (b.m.e.+A+B). In the last two columns we have added either the Jastrow (C) or UCOM (D) short 
range correlations. The value gx = 1-25 was used in these calculations. 

Table Hi] shows that the Jastrow method produces a much larger reduction in the magnitude of M^ ^ than the 
UCOM. The UCOM (D) with higher-order term (A) and finite nucleon size (B) corrections included seems to reduce 
the mag nitude of by a rough factor of 2/3 from its bare value. A similar scaling factor was also present for the 

OvP/3 results of 76 Ge and 82 Se in At the same time one obtains a reduction factor of 1/2 for the corresponding 
Jastrow (C) results. 

We visualize the <7 PP dependence of in Fig.[TJ The calculations included the higher-order term (A), finite-size 

(B) and UCOM (D) corrections. The calculations are shown for the interval g pp < 1.1, the upper limit lying near the 
breaking point of the pnQRPA for all shown nuclear systems. 

In Fig. [2] we have plotted the multipole decomposition of the total (V/?/3 matrix element M^ u ^ for all calculated 
decays. The upper (lower) end of each bar corresponds to the lower (upper) end of the g pp interval of Table H For 
the 136 Xe decay we can only give upper limits since only the lower limit of g pp is known as indicated in the last 
line of Table [B From the Fig. H one can see that the widest spread appears in the bar corresponding to the 1 + 



TABLE II: Total matrix element M' "' of J3|) computed by correcting the bare matrix element (b.m.e) for the higher-order terms 
of the nucleonic weak current (A), for the finite nucleon size (B), and for either the Jastrow (C) or UCOM (D) correlations. 
The used value of g pp is indicated in the first column. 



Nucleus gpp b.m.e. +A +A+B +A+B+C +A+B+D 

^Zr 1.085 -5.308 -4.814 -3.736 -2.454 -3.521 

100 Mo 1.08 -6.126 -5.571 -4.358 -2.914 -4.113 

116 Cd 0.99 -5.726 -5.172 -4.263 -3.169 -4.076 

128 Te 0.905 -7.349 -6.673 -5.260 -3.563 -4.979 

130 Te 0.87 -6.626 -6.021 -4.777 -3.285 -4.530 

136 Xe 0.74 -4.715 -4.269 -3.478 -2.537 -3.317 




FIG. 1: (Color online) Calculated values of M*- "' for the indicated 0vf3f3 decays as functions of g pp . The UCOM (D) was used 
with (?a = 1-25 including all the other (+A+B) corrections. 



contribution. Furthermore, the g pp interval extracted from the 2z//3/3 data confines the 1 + contribution in a striking 
way: for 96 Zr and 100 Mo the 1 + contribution is of opposite sign to the other contributions. This interference with the 
rest of the contributions reduces the magnitude of M^ ^ for these two decays. Another notable feature is that the 1~ 
contribution is always the leading one, the 2~ contribution being usually of comparable size. This pattern is different 
from the one of 76 Ge and 82 Se decays where the 2~ contribution was the dominant one [26]. All these observations 
are in qualitative agreement with the results of [13] ■ 

We have collected the obtained g pp limits of Table U and their corresponding 0^/3/3 matrix elements of ([3]) in 

Table IIIII It is worth pointing out that the magnitudes of the tensor matrix element M^ ^ are quite small and are 
not indicated in the table. As discussed earlier, the <7 PP limits arise from both the 2v/3[3 and beta decay data. For 
136 Xe only experimental lower limit of the 2v(3/3 half-life is available, yielding only a lower limit for <7 PP . For 128 Te the 
available beta decay data give g pp = 0.86 and for 116 Cd they give g pp = 0.82. In the last column of Table Hill we also 
tabulate our predicted half-life limits in units of yY/({m u ) [eV]) 2 . The beta-decay data yields a different g pp interval 
than the data on 2v(ij3 decays for 116 Cd and 128 Te. This discrepancy is especially striking for 116 Cd. Implications of 
this discrepancy and its cure are still open questions (2l| . 

Our results for Qvf3(3 nuclear matrix elements disagree with those of (2(J and (28|. In fact, just recently [36j the 
Tubingen-Caltec collaboration has corrected their results for a coding error in their computer program. The results of 
the Erratum [36j agree nicely with our results for the Jastrow corrected nuclear matrix elements. This means that one 
can safely say that the Jastrow short-range correlations reduce the values of matrix elements some 25%-40%. On the 
other hand, for the UCOM we obtain only a 4%-16% reduction. For this reason our present, UCOM corrected matrix 
elements are larger than the Jastrow corrected ones. Such differences give rise to big differences in the predicted Ovflfl 
half-lives for a given value of the effective neutrino mass (m u ). This invariably affects the sensitivity estimates for the 
presently running and planned double beta experiments. 
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FIG. 2: Multipole decomposition of M' "' for the calculated Of/3/3 decays. 



TABLE III: Calculated OuP/3 nuclear matrix elements, the used g pp and g\ values and the resulting half-lives. The UCOM and 
other corrections are included. The half-lives are expressed in units of yr/((m„)[eV]) 2 . 



Nucleus 


9pp 


9A 


Mp "' 


M GT 


M (tv) 


l l/2 


9B Zr 


1.06 


1.00 


1.350 


-2.969 


-4.319 


6.1 x 10 23 




1.11 


1.25 


1.261 


-2.315 


-3.117 


4.7 x 10 23 


luu Mo 


1.07 


1.00 


1.583 


-3.266 


-4.849 


6.2 x 10 23 




1.09 


1.25 


1.543 


-2.950 


-3.931 


3.8 x 10 23 


iib Cd 


0.82 (/T decay) 


1.25 


1.427 


-4.021 


-4.928 


2.3 x 10 23 




0.97 


1.00 


1.310 


-3.372 


-4.682 


6.3 x 10 23 




1.01 


1.25 


1.275 


-3.124 


-3.935 


3.6 x 10 23 




0.86 (/3~ decay) 


1.25 


1.939 


-4.276 


-5.509 


5.2 x 10 24 




0.89 


1.00 


1.866 


-3.975 


-5.841 


1.1 x 10 25 




0.92 


1.25 


1.792 


-3.650 


-4.790 


6.9 x 10 24 


130 Te 


0.84 


1.00 


1.699 


-3.743 


-5.442 


5.3 x 10 23 




0.90 


1.25 


1.575 


-3.219 


-4.221 


3.5 x 10 23 


13b Xe 


0.74 


1.00 


1.104 


-2.615 


-3.719 


1.1 x 10 24 



V. SUMMARY AND CONCLUSIONS 



We have calculated the Ov/3/3 nuclear matrix elements for the decays of 96 Zr, 100 Mo, 116 Cd, 128 Te, 130 Te and 
136 Xe by using the proton-neutron quasiparticle random-phase approximation with realistic two-body interactions in 
realistic single-particle spaces. We have corrected the bare matrix elements for higher-order terms of the nucleonic 
weak currents, for the nucleon's finite-size and for the nucleon-nucleon short-range correlations. The short-range 
correlations have been included by using the unitary correlation operator formalism. This method is superior to the 
rudimentary Jastrow procedure, traditionally adopted for the OV/3/3 calculations. 

The UCOM reduces the magnitudes of the matrix elements less than the Jastrow procedure. This leads to larger 
matrix elements and shorter 0vP(3 half-lives as compared to some recent calculations quoted in the literature. This 
has a notable influence on the estimated neutrino-mass sensitivities of the presently running and future double beta 
experiments. 
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